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Abstract 

A  linear  slender-body  theory  is  pre¬ 
sented  for  the  zero-speed  motions  of  a 
ship  in  regular  free-surface  waves, 
valid  for  an  arbitrary  angle  of  inci¬ 
dence  and  all  wavelengths  of  practical 
interest.  The  velocity  potential  of 
the  "radiation"  problem  is  obtained  by 
the  superposition  of  a  particular  solu¬ 
tion,  identical  to  that  of  the  short- 
wavelength  strip  theory,  and  a  homo¬ 
geneous  component  that  accounts  for  the 
longitudinal  flow  interactions.  A 
similar  representation  is  used  for  the 
"diffraction"  problem  where  the  particu¬ 
lar  solution  is  equal  and  opposite  to 
the  incident  wave  velocity  potential. 

Computations  are  presented  for  the 
vertical  hydrodynamic  force  distribution 
and  the  heave  and  pitch  added-mass  and 
damping  coefficients,  exciting  forces 
and  motions  for  a  Series  60  hull.  Com¬ 
parisons  are  made  with  strip  theory,  a 
three-dimensional  theory  and  experi¬ 
ments  . 


I,  Introduction 

In  the  past  twenty  years  slender- 
body  theory  techniques  have  found  broad 
applications  in  ship  hydrodynamics. 

The  geometrical  slenderness  of  the  ship 
hull,  however,  is  insufficient  to  jus¬ 
tify  an  asymptotic  theory.  The  length- 
scales  introduced  by  the  wavelength  and/ 
or  the  ship  forward  speed  need  also  to 
be  considered. 

Tne  wavelengths  encountered  in  an 
ambient  seaway  are  comparable  to  both 
the  transverse  and  the  longitudinal 
dimensions  of  conventional  ships.  It 
is  therefore  desirable  that  a  theory 
for  the  motions  of  a  ship  advancing  in 
waves  would  embrace  all  wavelengths 
encountered  in  practice.  Newman  (1978) 
developed  the  theoretical  framework  of 
such  a  slender-body  theory  where  the 
ship  slenderness  is  sufficient  for  the 
asymptotic  analysis  involved.  The 
theory,  hereafter  referred  to  as  "uni¬ 
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fied  theory",  was  applied  to  the  radia¬ 
tion  problem  bv  Newman  and  Sclavounos 
(1980)  and  was  extended  to  the  diffrac¬ 
tion  problem  by  Sclavounos  (1981). 

In  the  present  paper  the  radiation 
and  diffraction  problems  are  combined 
to  predict  the  zero-speed  heave  and 
pitch  motions  of  a  shin  in  waves.  The 
theoretical  derivation  is  presented  in 
Section  2  including  the  radiation  prob¬ 
lem  as  derived  by  Newman  (1978)  for 
purposes  of  comparison  with  the  diffrac¬ 
tion  problem.  The  short- wavelength 
approximations  of  the  unified  theory 
are  derived  in  Section  3  where  existing 
theories  for  the  radiation  and  diffrac¬ 
tion  problems  are  recovered.  An  effi¬ 
cient  numerical  technique  for  the  solu¬ 
tion  of  two-dimensional  free-surface 
problems  is  described  in  Section  4.  A 
hybrid  integral  representation  valid  in 
a  domain  bounded  by  the  body  section, 
the  free  surface  and  a  circular  boundary 
lying  in  the  fluid  domain  is  matched  to 
a  multipole  representation  valid  outside 
the  matching  circle.  In  Section  5  numer¬ 
ical  computations  are  presented  for  the 
hydrodynamic  force  distribution  and  the 
motions  for  a  Series  60  (Cgra0.7)  hull 
and  comparisons  are  made  with  existina 
theories  and  experimental  data. 

In  what  follows, a  physical  descrip¬ 
tion  of  the  method  of  solution  is 
attempted.  The  ship  is  assumed  to  be 
slender,  the  fluid  motion  incompressible 
and  irrotational  and  the  ship  motion 
amplitude  small  enough  to  validate  the 
linear  decomposition  of  the  hydrodynamic 
disturbance  in  the  radiation  and  the 
diffraction  problems. 

For  the  radiation  problem,  in  an 
"inner"  region  close  to  the  ship  hull, 
the  longitudinal  flow  gradients  are 
small  compared  to  the  transverse  ones. 

The  field  eouation  can  therefore  be 
reduced  to  the  two-dimensional  Laplace 
equation, supplemented  by  a  two-dimen¬ 
sional  bodv  boundary  condition  and  the 
zero-speed  free-surface  condition.  For 
the  diffraction  problem  the  slenderness 
assumptions  are  invoked  after  the 
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longitudinal  component  of  the  incident 
wave  is  factored  out.  The  two-dimen¬ 
sional  modified  Helmholtz  equation  turns 
out  to  be  the  relevant  inner -field  equa¬ 
tion  with  a  no  flux  body-boundary  condi¬ 
tion  and  the  same  free-surface  condition 
as  in  the  radiation  problem. 

In  an  "outer"  region/  located  at 
radial  distances  greater  or  equal  to  the 
ship  length/  the  flow  gradients  in  all 
directions  are  of  comparable  magnitude. 
The  three-dimensional  Laplace  equation 
is  satisfied  in  the  fluid  domain  for 
both  problems/  subject  to  the  zero- 
speed  free-surface  condition  and  a  con¬ 
dition  of  outgoing  waves  at  infinity. 

The  method  of  matched  asymptotic  expan¬ 
sions  is  used  to  enforce  the  compatibi¬ 
lity  of  the  inner  and  outer  solutions 
in  an  intermediate  "overlap"  region. 

The  velocity  potential  of  the  inner 
radiation  problem  involves  a  particular 
solution  identical  to  that  of  the  short- 
wavelength  strip  theory  and  is  supple¬ 
mented  by  a  homogeneous  solution  that 
represents  standing  waves  at  "infinity". 
The  corresponding  particular  solution 
of  the  diffraction  problem  is  equal  and 
opposite  to  the  incident  wave  potential 
with  the  longitudinal  wave  component 
factored  out.  The  homogeneous  part 
again  satisfies  standing  waves  at 
"infinity"  and  is  regular  for  all  head¬ 
ings. 

The  homogeneous  solutions  of  both 
the  radiation  and  the  diffraction 
problems  are  multiplied  by  "interaction" 
coefficients,  functions  of  the  longitu¬ 
dinal  coordinate,  that  are  determined 
after  the  matching  with  the  correspond¬ 
ing  outer  solutions.  Two  integral 
equations  are  solved,  in  that  context, 
for  the  unknown  three-dimensional  source 
distributions  along  the  ship  axis  that 
describe  the  outer  solutions  of  both 
problems. 

In  the  short-wavelength  limit,  the 
interaction  coefficient  of  the  radiation 
problem  vanishes  and  the  classical  strip- 
theory  solution  is  recovered  (  see 
Ogilvie  and  Tuck(1969)).  At  the  same  limit, 
the  solutions  of  Choo  (1975)  and  Troesh 
(1976)  for  the  diffraction  problem  are 
recovered  for  oblique  waves.  The  sin¬ 
gularity  that  is  present  in  these 
theories  for  head  waves  is  due  to  a  non¬ 
uniformity  that  is  present  when  the 
short-wavelength  approximation  preceeds 
the  transition  from  oblique  to  head 
incidence.  For  head  waves  the  Maruo 
and  Sasaki  (1974)  theory  is  recovered, 
apart  from  a  minor  analytical  differ¬ 
ence. 


2.  The  Boundary  Value  Problem 


We  introduce  a  Cartesian  coordinate 
system  fixed  in  space  with  the  free 
surface  taken  at  z=0,  the  ship  center- 
plane  at  y=0  and  the  positive  axis  point¬ 
ing  towards  the  bow.  The  ship  is 
assumed  to  perform  small  oscillatory 
motions  around  its  undisturbed  floating 
position,  all  oscillatory  quantities  are 
expressed  in  complex  form  and  the  factor 
elljt  is  understood  hereafter.  Only 
the  symmetric  heave  and  pitch  motions 
are  considered,  since,  to  the  same 
degree  of  accuracy,  the  longitudinal 
interactions  are  absent  from  the  anti¬ 
symmetric  modes  for  which  strip  theory 
is  valid  for  all  wavelengths.  More 
details  for  the  analysis  that  follows 
can  be  found  in  Newman  (1978)  and 
Sclavounos  (1981). 

For  water  of  infinite  depth  the 
incident  wave  potential  is  given  by 


*  -3£  exp[vz-iv (xcosfi+ysinp) )  ,  (2.1) 

where  A  is  the  wave  amplitude,  v=(dj/ g  is 
the  wavenumber  and  B  the  angle  of  inci¬ 
dence,  with  6*180°  for  head  incidence. 

Let^.,  j*3,5,7  be  the  complex 
velocityvotentials  associated  with  the 
unit  ampli -ude  heave  and  pitch  motions 
and  the  symmetric  diffraction  problem 
respectively.  Under  the  assumptions 
stated  in  the  Introduction,  the  three- 
dimensional  Laplace  equation  is  satis¬ 
fied  in  the  fluid  domain 


>j  *  °' 


subject  to  the  linear  free-surface  condi¬ 
tion 


^jz”Vfj 


on  z  *  0  , 


and  the  body  boundary  condition  applied 
at  the  mean  position  of  the  ship  vetced 
surface  S’ 


,■*  vz-ivxcosP 

(n.V)  (---  e 


,j»3,5 

cos  (vysinP) ) ,  j=7 


.  (2.4) 
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where  n  =  (n.,n.,nJ  is  the  unit  normal 
vector  pointing  out  of  the  fluid  domain 
and  n5  =  -xn,+zn..  The  well-posedness 
of  the  previous  Boundary-value  problems 
requires  that  the  velocity  potentials 
represent  outgoing  waves  in  the  far 
eld  with  the  fluid  velocity  vanishing 
as  z-*— «. 

In  the  next  two  subsections  the 
ship  slenderness  is  invoked  to  approxi¬ 
mate  the  flow  equations  in  the  outer  and 
inner  regions  respectively. 


2A.  Outer  Region 


For  small  vr,  vQr  and  kr,  G  can  be 
approximated  in  the  form 

G*(y,z;k-v0)  =  Re(G* (y,z;- JQ) } 

-  yjj-d+vz) F* (k-vQ)  (2.8) 

with  an  error  factor  l+0(v2r2,v  2r2,k2r2) 
and  ° 

F*  (k)  =  £n{ |vQ|/|k | )  +  (l-vo*/v2  )-5i 


At  radial  distances  comparable  to  or 
greater  than  the  ship  length,  the  flow 
is  insensitive  to  the  ship  hull  geometry 
details  and  the  velocity  potential  can 
be  approximated  by  a  line  distribution 
of  three-dimensional  sources  along  the 
ship  centerline,  given  by 

iv  x 

*  e  <J>j(x,y,z) 


-  J  e^  qj(^)  G(x-C,y,r)dC,  (2.5) 
L 


where  v  is  the  characteristic  wavenum¬ 
ber  of  the  longitudinal  distribution  of 
the  hull  normal  velocity  (vo«0  for  the 
heave  and  pitch  motions  and  v  »  -  vcosB 
for  the  diffraction  problem) ,  qj (x)  is 
the  slowly  varying  part  of  ths  source 
strength  distribution  and  G(x-C,y,z)  is 
the  velocity  potential  due  to  a  unit 
pulsating  source  located  on  the  ship 
axis  at  x=f. .  In  order  to  match  the 
outer  solution  (2.5)  to  the  correspond¬ 
ing  expression  of  the  inner  problem  we 
need  to  expand  (2.5)  for  small  vr, 
where  r  =  (y2+z2)*5.  It  is  convenient 
to  take  the  Fourier  transform  of  both 
sides  of  (2.5)  to  obtain, 


•J'j*  (y,z;k)  =q//k)  G* (y, z;k-vQ) ,  (2.6) 


•  cosh-1  ( | v/  vo|)  -  |  l-k2/v2  j-1* 


ni+cosh-1 (v/ |k | ) 
-tt+  cos-1  (v/|k  | ) 


where  the  upper  or  lower  term  in  brack¬ 
ets  is  applicable  according  as  v/|k|£l 
and  |v  |<_  v  for  all  cases  of  interest. 


Expression  (2.8)  is  a  generalization 
of  the  corresponding  approximation 
derived  bv  Newman  {1978,  eq.  4.12) 
for  the  radiation  problem.  As  v  *0, 
the  combination  of  the  first  two°terms 
in  (2.9)  reduces  to  £n(2v/|k|). 
Furthermore,  G  iy,z>0)  *  R_  (y,z)  is 
the  two-dimensional  source^potential 
that  satisfies  the  Laplace  equation, 
subject  to  the  linear  free-surface 
condition  and  a  radiation  condition 
of  outgoing  waves  qs  v|y|  •»  .  Given 
that  ImtRjp)  *  ^  e'*cos(vy)  and  the 

fact  that  evzcos(vy)  ■  (1+vz)  (J.+0  (v2r2 )  1 
in  the  regime  where  (2.8)  is  valid,  it 
follows  that 


G*  (y,z  jk)  *  R2D(y,i-  )-  y^(l+vz)  f*  (k)  , 

(2.10) 


where 


* 

G 


(y,z;k) 


(«»>  , 

^  eikx  G(x,y,z)dx 


z  (u‘+k  ) 
e  cost 

(u2+  k‘ 


where 

f*  (k)  -  fn(2v/|k|)  +  ni  -  1 1-k  ’  /v2 |  . 

ni+cosh"1 (v/ | k | fj  ^  j 

-*  +  cos-*  (v/|k  | )  ! 


(2.7) 


which  is  identical  to  equation  (4.12) 
of  Newman  (1878). 


dinate  stretching  suggests  that 


For  the  diffraction  problem, vQ  = 

-vcosB,  with  G*(y,z;-v  )  =  D2d  being 
the  corresponding  two-dimensional 
source  potential  that  satisfies  the 
Helmholtz  equation  (V|D-v02)D2D  =  ® 

and  the  same  free-surface  and  radiation 
conditions  as  R~_.  It  is  essential  to 

1  vz 

point  out  that  Im(D2D)=jcsci3e  cos(vysinB) 

is  singular  for  ]cosB|=l  and,  as  opposed 
to  the  radiation  problem,  it  is  not  present. 
Thus,  the  inner  expansion  of  the  outer 
solution  for  the  diffraction  problem 
is  given  by  (2.6)  and  (2.8)-(2.9)  and 
is  regular  for  all  headings. 

In  summary,  the  inner  expansion  of 
the  outer  solution  transformed  back  in 
the  physical  x-space  takes  the  form 


i|<j(x,y,z)  *  q^(x) 


Re(D2D) 


^j-d+vz)  /.  (qj) , 


(2.12) 


where  the  linear  operator  JL  is  obtained 
by  taking  the  inverse  Fourier  transform 
of  (2.6).  It  follows  that 


^(q.)  =  q.(x)  [v  +  {  h(6)}] 


-ix0(x-n  i 


{ ^sgn (x-C)fn(2v|x-E | ) 


TS-  0(1>'  =  0<£_1) 


(2.16) 


when  operating  on  defined  in  (2.5). 

Applying  (2.16)  in  (2.2)-(2.4)  and 
using  the  definition  of  <Jj., we  obtain 
the  two-dimensional  modified  Helmholtz 
eauation 


\ i> .  +  <J> .  —  v  2tl> .  =  0  , 

iyy  o 


(2.17) 


subject  to  the  free-surface  condition 


<|i.  -  .  =0  on  z=0  , 

1 z  j 


and  the  body  boundary  condition 


(2.18) 


W  *j’ 


iu  N.  ,  j=3, 5 

J  r> 

-3-  -^{evzcos(vysinB)  },  j=7 


(2.19) 


where  to  leading  order  in  e,n*  (0,N.N,) , 
with  H5*  -  xN3. 

The  general  solution  of  (2.17)- (2. 19) 
can  be  obtained  in  the  form 


Vj  =  ♦jp  +  Cj(X)  *jH» 


(2.20) 


’((Sr+  ivo)qj <f,)  ■  T  K(  '»(x-C)]qj(C)»» 

(2.13) 

where  y-Q.Sl...  is  the  Euler  constant, 
h (6)  =cscBcosh  1  ( | sec8 1 ) -£n(2|sec8 1)  (2.14) 


respectively.  The  "interaction" 
coefficient  C • (x)  is  presently  an 
arbitrary  function  of  x  that  will  be 
determined  from  the  matching  with  the 
outer  solution.  No  radiation  condition 
needs  to  be  satisfied  in  the  inner 
region  since  the  outer  solution  already 
represents  outgoing  waves  at  infinity. 


K(x)=YQ(x)+2i  J0(|x|)+  «Q ( | x j )  (2.15) 


with  the  upper  or  lower  terms  in  brack¬ 
ets  applying  for  the  radiation  (j=3,5 
and  v<j=0)  and  the  diffraction  (j=7  and 
v  =  -vcosB)  problems  respectively. 


Inner  Region 


At  transverse  distances  of  the  order 
of  the  ship  beam,  the  relative  order  of 
the  flow  gradients  in  the  longitudinal 
and  transverse  directions  are  dictated 
by  the  respective  order  of  the  geometry 
gradients.  Assuming  that  both  the  y 
and  z  coordinates  are  of  0(0,  a  coor- 


Following  Newman  (1978),  a  particular 
solution  of  the  radiation  problem 
(j=3,5,  v  *0)  is  the  short-wavelenqth 
strip  theory  potential  4>j  that  satisfies 
outgoing  waves  as  vjy |  The  homo¬ 
geneous  solution  can  be  physically 
regarded  as  the  interaction  of  two  inci¬ 
dent  symmetric  waves  with  the  ship 
section.  The  pure  imaginary  form  of 
the  body  boundary  condition  suqgests 
that  4>j  +  is  a  homooeneous  solution, 

where  the  overbar  stands  for  the  com¬ 
plex  conjugate  of  the  quantitv  involved. 
Thus, 

=  *.+Cj(x)  (fy  J.),  j  =  3,5  (2.21) 


the  ship  section,  and  is  obtained  in 
the  form 


with  In  the  overlap  region,  at 

transverse  distances  large  compared  to 
the  ship  beam  but  still  small  compared 
to  its  length  as  e+0,  can  be  ex¬ 
pressed  by  their  effective  source 
strengths,  in  the  form 


♦j  =  °j  R2D(y'2)  ' 


(2.22) 


where  R2D  is  the  two-dimensional  source 
potential  that  satisfies  outgoing  waves 
as  v|y|-><».  We  may  thus  approximate  A  . 
in  the  overlap  region  as  follows,  '  3 

h  «  {aj+Cj<x)<oj+ai)}R2D 


-  2i  Cj (x)Oj  Im(R2D). 


(2.23) 


If  we  now  set,  Im(R2D)*i(l+uz)  (l+0(v2r2)}, 
the  outer  expansionist  the  inner  solu¬ 
tion  becomes 

1 } j  *  {Oj+C^xHOj+o.))  R2D 


#7H  “  eVZ  cos(vysinB)+<(iD(y,z),  (2.29) 


where  \pD  is  constructed  by  using  the 
two-dimensional  source  potential  Re(D2D) 
that  represents  standing  waves  as 
v|y|-»«  and  is  regular  for  all  headings. 
Thus, 


^7  s  ~  u  cos (vysinS) (1-C7 (x) 1 


+  C7(xHD(y,2). 


(2.30) 


In  the  overlap  region,  can  be 
approximated  in  the  form 

<J>7  *  -  ^  ll-C7(x) )  (l+vz) 


+  C7(x)o?(x)  Re(D2D) 


(2.31) 


-  i  Cj (x) o j (l+vz). 


(2.24) 


Comparing  (2.12)  to  (2.24)  the  match¬ 
ing  conditions  for  the  radiation  prob¬ 
lem  are. 


Oj  +  Cj (x) (Oj+Oj)  -  qj 

i  Cj  ^<qj>* 


(2.26) 


After  eliminating  Cj  from  (2.25)-(2.26) , 
the  outer  source  strength  q.  is  deter¬ 
mined  from  the  integral  equation 

qj(x)  “  5Fi (oj/°j+1)^(qj) 


*  Oj  (x)  j*3, 5  , 


(2.27) 


where  *(qj)  is  given  by  (2.13)  with  vQ*0» 

A  particular  solution  for  the  dif¬ 
fraction  problem  follows  easily  from  the 
body  boundary  condition  (2.19)  and  is 
equal  and  opposite  to  the  symmetric 
part  of  the  incident  wave  potential 


The  matching  requirements  again 
follow  easily  from  (2.12)  and  (2.31), 


C7  °7  *  q7 


(2.25)  &U-C7)  - 


(2.32) 


(2.33) 


The  elimination  of  C7  from  (2. 32)- (2. 33) 
results  in  the  following  integral 
equation  for  q7, 

<*7(x)  -7£87  *£(<J7)m  °7  '  (2’34> 


where  ^(q7)  is  < 
with  v0»-vcos(J. 


determined  by  (2.13) 


The  solution  of  equations  (2.27)  and 
(2.34)  determines  the  outer  source 
strength  q.,  and  the  complete  inner 
solutions  follow  in  the  form 


V'i 


(♦>♦<>  j*3,5  (2.35) 


for  the  radiation  problem,  and 


<^7p  =  -  ^  evz  cos(vysinB).  (2.28)  ^7  “  eVZ  cos(vysinB) (q7/a7  -  1) 


The  head  sea  limit  of  is  the  leading 
order  solution  in  the  snort-wavelength 
theory  of  Faltinsen  (1971). 

The  homogeneous  solution  again  re¬ 
sults  from  the  interaction  of  two  waves 
of  unit  amplitude  and  equal  phase  in¬ 
cident  from  opposite  directions  upon 


+  *D  V°7 


for  the  diffraction  problem. 


(2.36) 


3.  Short-Wavelength  Approximations 

For  wavelengths  comparable  to  the 
ship  transverse  dimensions  |vlr=0(e-i)  as 
e-*0),  it  is  possible  to  approximate  the 
linear  operator  g£(a)  defined  in  (2.13) 
using  the  techniques  of  asymptotic 
analysis.  For  the  diffraction  problem 
and  -l<cosB£0,it  follows  that 

.  ,  .  V2  -iv(£-x)(l+cosS) 


<5  -  x) 


7)  =  j(l-i)  (nv*5)  Jq7(£) 


where  1+vz  =  ev  cos(vysinB) (l+0(v2r2) J 
in  the  overlap  region.  The  two-dimen¬ 
sional  source  potential  D2d  satifies 
outgoing  waves  as  v|y|-*<*>  is  singular 
for  head  waves  and  corresponds  to  the 
Green  function  used  in  the  short-wave¬ 
length  strip  theory  derived  by  Choo 
(1975)  and  Troesh  (1976). 

In  the  head-sea  case  (6=180°),  (3.1) 
reduces  to  (originally  due  to  'h.itinsen 
(1971))  L/2 


-7?i (cscB-  ^|sec  || )q? (x)+0 (v'S.  (3.1)  =  J  <l-i) <»v)*  j 


The  details  of  the  derivation  can  be 
found  in  Sclavounos  (1981).  The  corre¬ 
sponding  approximation  for  the  radiation 
problem  follows  from  (2.13)  and  (3.1) 
with  cos6=0,  in  the  form 

»  1  uY2  e-iv(*-x) 

iC(q.)=i(l-i)  (ttv)*  q.(£)- - - 

3  2  l  3  <e-xP 

+  ni  ^  qj(x)  +  O^**).  (3.2) 

Proceeding  with  the  radiation  problem, 
the  integral  in  (3.2)  can  be  further 
approximated  for  large  v  in  the  form 


J  i 


-iv(C-x) 


dC  =  -<^)Vli/4qj(x). 


The  substitution  of  (3.3)  in  (3,2)  can¬ 
cels  out  the  first  two  terms,  with  the 
remaining  part  being  of  0(v“’).  Thus 
for  short  wavelengths 

q  j  (x)  =  Oj(x)  +  0  ( v*”*5 )  ,  (3.4 

and  strip  theory  is  recovered. 

An  approximation  similar  to  (3.3) 
applied  to  (3.1)  gives 

^(q?)  =  -Hi  cscB  q7  (x)  +  0  (v-*5 ) .  (3.5 

The  approximation  (3.5)  is  singular 
for  head  waves  due  to  a  nonuniformity 
that  is  present  in  the  asymptotic 
approximation  of  (3.1)  when  l+cos8*0. 
Substituting  (3.5)  in  (2.34)  we  obtain 
an  algebraic  equation  for  q7(x)  with  a 
solution 

q7  =  o?/(l  +  o7cscB) .  ( 3, ( 


Substituting  (2.33)  and  (2.32)  in 
(2.30)  and  using  (3.5),  we  may  write 
the  outer  expansion  of  the  inner  solu¬ 
tion  in  the  form 


and  the  Maruo  and  Sasaki  (1974)  theory 
is  obtained  apart  from  two  analytical 
differences;  the  *j  factor  is  missing 
in  their  theory  and  an  additional 
Hi  q7 (x) /2  term  appears  in  their  eaua- 
tion  which  corresponds  to  (3.8)  that 
should  not  be  present  in  a  consistent 
short-wavelength  approximation. 


4.  The  Two-Dimensional  Problem 

In  this  section  a  numerical  scheme 
for  the  solution  of  two-dimensional 
free-surface  wave  problems  in  water  of 
infinite  depth  is  described. 

The  application  of  spectral  tech¬ 
niques  to  the  solution  of  free-surface 
boundary-value  problems  is  in  general 
restricted  to  simple  body  profiles. 
Ursexl  (1949)  obtained  the  solution 
of  the  heaving  problem  of  a  circular 
section  in  water  of  infinite  depth  by 
using  a  sequence  of  singularities 
located  at  the  origin  of  the  axes. 

The  extention,  however,  of  this  method 
to  more  general  body  sections  has  not 
been  yet  formally  established. 

An  efficient  numerical  technique 
to  treat  arbitrary  body  and  bottom 
geometries  has  been  developed  by  Yeung 
(Bai  and  Yeung  (1974)].  In  an  inner 
region  of  changing  topography,  a  hybrid 
integral  equation  is  obtained  by  apply¬ 
ing  Green's  theorem  and  usinq  the  funda¬ 
mental  logarithmic  singularity  as  the 
relevant  Green  function.  This  repre¬ 
sentation  is  then  matched  to  an  eigen¬ 
function  expansion  in  an  outer  region 
of  constant  depth.  Liapis  and  Faltinsen 
(1980)  extended  this  method  to  the  case 
where  the  Helmholtz  eouation  is  satis¬ 
fied  in  the  fluid  domain.  A  compre¬ 
hensive  survey  of  existing  techniques 
for  the  solution  of  two-dimensional 
free-surface  wave  problems  is  given  by 
Wehausen  (1974). 


=  j  i  cscB  q_(x)  (1+vz) +q? (x) Re (Djjj) 
=  q?(x)  D2d  ,  (3. 
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where  K  (x)  is  the  modified  Bessel 
function  of  order  n  defined  in  Abramowitz 
and  Stegun  (1964).  For  £=0,  W  =R»n 
with  R2D  defined  in  Section  2,andz 
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Applying  Green's  theorem  in  region 
II  (ABCDA) ,  using  the  symmetry  of  the 
flow  and  the  free  surface  condition 
(4.2),  we  obtain 


The  present  approach  is  a  coupling 
of  a  multipole  expansion  and  a  hybrid 
integral  representation,  valid  respec¬ 
tively  in  the  domains  I  and  II  of  Figure 
1,  separated  by  a  matching  circle  SR. 

The  water  depth  is  assumed  infinite  and 
the  body  section  symmetric  with  respect 
to  the  y«0  axis.  The  Helmholtz  equation 
is  in  general  satisfied  in  the  fluid 
domain 

Va4>  -  £2 $  *  0,  (4.1) 

subject  to  the  free-surface  condition 


4>  -  v  <J>  «  0  on  z*0,  (4.2) 

z 

to  a  general  body-boundary  condition 

(n  •  V ) <j>  »  V(P)  on  Sn  ,  (4.3) 

P  B 


-l$n(Q)  +  [  dstPH11  (P)  (n  -V  )F(P,Q) 

£  I  P  P 

SBUSR 

+  J  ds(P)^II(P)l(np-Vp)-v)F(P,Q) 

SF 

-  J  ds(P)F(P,Q)  (np‘Vp)4,1I(p) 

SR 

«  |  dS(P)  F(P.Q)  V(P),  (4.7) 

SB 

where  np  is  the  unit  normal  vector 
pointing  out  of  domain  II  and 


and  to  a  condition  of  standing  or  out¬ 
going  waves  as  v|y|-*®.  For  all  cases  of 
interest  0<^  t  v . 

Assuming  a  symmetric  flow  with  re¬ 
spect  to  the  y*0  axis,  we  may  write 
the  velocity  potential  in  region  I  in 
the  form  (Urnell  (1949,  1968)] 


F(P,Q)  -  G (P,Q  +  )  +  G(P,<T),  (4.8) 

with  Q+=  Q  lying  on  SgUSpUSR  and  Q- 

being  the  symmetric  of  Q  with  respect 
to  the  y»0  axis;  G(P,Q)  is  the  velocity 
potential  at  P  due  to  a  fundamental 
unit  source  located  at  Q,  given  by 


aQ  Wo(r,0) 


m*l 


am  Vr'0)' 


(4.4) 


where  aQ  is  the  effective  source 
strength,  WQ  is  the  relevant  wave 
source  potential  due  to  a  unit  source 
located  at  the  origin, and  W„(r,9)  are 
wave-free  multipoles.  For  4*  v|cosB|>0, 
W0»Re(D2I))  with  D2d  defined  in  Section 
2 ,  and 


G(P,Q  ) 


(’nr*°g(rpq)  * 
^?rKo(trpQ)  '  *>0 


(4.9) 


where  K  (x)  is  the  modified  Bessel 
function  of  zeroth  order. 


A  grid  is  selected  along  the  contours 
S„,  S  and  SR  with  the  velocity  potential 
assumed  constant  along  the  segments 
joining  two  consecutive  grid  points. 

Let  N  ,  N_  and  Nr  be  the  number  of  seg- 
mentsBalohg  S_,  SF  and  S_  respectively. 
Taking  the  segment  lengths  to  be  of 
equal  length  along  the  matching  circle 
and  keeping  N„-l  terms  in  the  series 
(4.4),  the  velocity  potential  and 
its  r-derivative  41  at  the  segment 
midpoints  can  be  written  in  a  matrix 
form 


I  ** 
r  *  A  a 


$r  -  B  a  , 


(4.10) 

(4.11) 


where  the  column  vectors  9  ,  a  and  the 
matrixes  A  and  B  are  defined  as  follows 


where  NT  =  N„+NF+N  ,  V.  is  the  body 
boundary  normal' velocity ,  6. .  the 


Kronecker  delta  and 


£>..  j—  1, . •  ,Njj  and 3“N— ■ +1, . . 

C.j-  13  -B  7  T  (4>17) 

Dij-^ij  3=^+1,..,  Np 

with  i=l,..,N  .  The  influence  matrixes 
S. D. .  i,j  *  1,...,N  are  defined  as 
f&Ilowi3  T 


Sij  =  {  ds (P)F(P»Qi>  (4.18) 

Sj 

Dij  =  |  ds(p)  (np-Vp)F(P,Qi),  (4.19) 


41  =  {4I(R,e1),...,4>I(R,eN  ) J1 

R 


a  =  {aQ,  • . • ,  aR 

R 


(4.12) 

(4.13) 


A  *  A.j=Wj-_1(R,ei),B=^  | 

i,j  *  1,...,NR,  (4.14) 

where  (R,6,),  i=l,..,,N_  are  the  polar 
coordinates  of  the  segment  midpoints 
along  the  discretized  circle  SR,  with 
0^  decreasing  to  zero  as  i  increases. 

If  we  left-multiply  both  sides  of 
equation  (4.10)  by  BA-1,  we  obtain  a 
relation  between  the  velocity  potential 
4>1  and  the  normal  velocity  ,  in  the 
form 


where, F(P,Q)  is  defined  in  (4.8)- (4.9) , 
the  point  Q.  is  located  at  the  midpoint 
of  the  i-th  segment  and  the  inteqrations 
in  (4,18)-(4,19)  are  performed  over  the 
j-t.i  segment  S^. 

The  matching  conditions  of  continuous 
pressure  and  normal  velocity  across  the 
matching  circle  S_  can  be  expressed  in 
the  form  K 


(4.14)  ^(R^.)  « 


4>  (R.6.)  *  (t11  . 
vrv  '  i  vr  , i+N 


(4.20) 


(4.21) 


where  i*l,...,NR  and  N=N„+NP.  Using 
(4 . 20) - (4 . 21)  and  (4.15) Bin  the  second 
sum  of  (4.16),  we  end  up  with  a  matrix 
equation  for  j*l,...,Ny 


*r  *  B  A  S  .  (4.15) 

(C-|  I-E}$n  =  S  V, 

(4.22) 

We  next  proceed  to  the  discretization 
of  equation  (4.7).  V.ie  velocity  poten¬ 

where , 

tial  at  the  midpoint  of  the  i-th  seg¬ 
ment  is  denoted  by  ,  where  the  in¬ 
dexing  starts  from  trie  point  A  of  Fig. 

411  *  {*JX)T  i-l,...,NT 

(4.23) 

1.  It  follows  that 

v  =  (vi,  0,..,o|  i“l, . . . ,Ng 

(4.24) 
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i,  j  *  1,...,  Nij.  (4.25) 
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E=Eij= 


k=N+lSik  m=i  k_N'm  n»,j-N, 
j=l i • • » ,N 


j  — N+l  i  ■  •  »  ,N(. 


(4.26) 


with  i=l,...,NT  and  A=Aij,  B=Bi ^ , 
i,j=l,...,NT  defined  in  (4.14). 

The  velocity  potential  $  along  the 
contours  SB,  Sp  and  SR  is  determined 

from  the  solution  of  the  system  of 
equations  (4.22).  The  effective  source 
strength  a  in  (4.4)  is  obtained  by 
left-multiplying  (4.10)  by  A“*. 

The  present  scheme  is  free  from 
irregular  frequencies.  Furthermore, 
only  one  evaluation  of  the  wave-source 
potential  WQ  is  necessary.  This  reduces 
substantially  the  total  computation  time 
required,  in  comparison  to  the  alterna¬ 
tive  approach  of  a  wave-source  distri¬ 
bution  on  the  body  boundary,  especially 
when  f >0. 

The  described  technique  has  been 
applied  for  f=v|cos6|>0  in  the  context 
of  the  diffraction  problem.  Selecting, 
the  matching  circle  radius  R  10%  greater 
than  the  body- section  maximum  radius  and 
15,  10  and  15  segments  on  SR,  Sp  and  SR 

respective^  proved  to  be  sufficient 
v  B 

for  0  <  -j-  <.  n ,  where  B  is  the  beam  of 
the  body  section.  The  same  estimates 
are  expected  to  hold  in  the  special 
case  of  t- 0.  For  this  paper,  the  two- 
dimensional  solutions  of  the  radiation 
problem  are  outained  by  using  a  com¬ 
puter  program  due  to  Yeung  (1975), 
based  on  his  hybrid  integral  equation 
technique. 


5.  Hydrodynamic  Forces  and  Motions 


where  is  the  complex  amplitude  of 
the  ship's  displacement  and 


Vi 


w  a.-  -  iu  b..  -  c.. 
ID  ID  ID 


(5.2) 


The  coefficients  a..,  b..,  c. .  are  real 
and  correspond  respectively  to  the  force 
components  in  phase  with  the  accelera¬ 
tion,  velocity  and  static  displacement 
of  the  ship.  The  added-mass  (a^)  and 

damping  (b. .)  coefficients  can  be 
derived  fr&m  the  inner  velocity  potential 
(2.21)  of  the  radiation  problem  by  inte¬ 
grating  the  linearized  pressure  obtained 
from  Bernoulli's  eauation  over  the  ship 
wetted  surface.  It  follows  that 


u^a .  .  -  iwb . . 
ID  ID 


—  ld)p 


|  ni$jds 


nX.  (x)  (41  j+41  j )ds, 


(5.3) 


where  all  the  Quantities  are  defined  in 
Section  2B. 

The  first  integral  in  (5.3)  corre¬ 
sponds  to  the  strip-theory  contribution. 
The  second  inteoral  is  the  correction 
due  to  the  three-dimensional  interaction 
effects  introduced  by  the  unified  theory. 
As  v+®,  it  follows  from  (2.25)  and  (3.4) 
that  C. (x)+0  and  the  strip-theory  is 
recoveted. 

Proceeding  to  the  diffraction  prob¬ 
lem,  the  resulting  complex  amplitude 
of  the  wave  exciting  force  is  simply 

Fi  =  A  X.,  i=3, 5 ,  (5.4) 


where  A  is  the  amplitude  of  the  incident 
wave.  Using  the  incident-wave  potential 
(2.1)  with  A=1  and  the  diffraction 
potential  (p the  unit-amplitude  ex¬ 
citing  force  X,  can  be  expressed  in  the 
form 


X;  =  -nop  I 
1  } 


+?7)  nids' 


(5.5) 


The  quantities  of  interest  for  the 
evaluation  of  the  ship  motions  are  the 
hydrodynamic  pressure  force  and  moment 
acting  on  the  ship  in  the  radiation  and 
the  diffraction  problems.  For  a  ship 
undergoing  a  steady-state  small  ampli¬ 
tude  oscillatory  heave  and  pitch  motion 
in  an  otherwise  calm  free  surface,  the 
resulting  complex  force  amplitude  can  be 
written  in  the  form 


An  alternative  expression  can  be 
obtained  by  using  the  Haskind  relations. 
Combining  the  body-boundary  condition 
(2.4) for  the  radiation  problem  and 
applying  Green's  theorem  in  (5.5),  we 
obtain 


Hi  =  )' 

j=3,5 


n  .  t .  , 

D  13 


1=3,5 


where  0» .  is  the  radiation  potential  due 
(5.1)  to  the  l-th  mode  of  oscillation.  With 

the  velocity  potentials*^  and ^  assumed 

known,  no  assumptions  are  involved 


• . 


V 


regarding  the  wavelength  or  the  body 
geometry  in  (5.5)  and  (5.6).  The 
advantage  however  in  using  (5.6)  is 
that  the  solution  of  the  diffraction 
problem  is  not  needed  for  the  evaluation 
of  the  exciting  force  and  moment  of  a 
body  in  waves. 

For  a  slender  ship,  the  inner 
diffraction  velocity  potential  defined 
in  (2.5)  and  (2.36)  is  used  together 
with  (2.1)  to  give 


Cj,  =  -iiopj  i 


-ivxcosB 


C?(x)  i(<7Hnids  (5.7) 


where  C -  and  <|>7H  are  defined  in  Section 
2B.  The  corresponding  expression  for 
(5.6)  is  obtained  by  substituting  the 
unif ied-theory  radiation  potential 
(2.21)  in  (5.6).  If  follows  that 

rr  afi 

xi  =  -Pjpun.fi  -  ♦i  -5-  )ds 


>11  ci 


(x) (♦,+?. ) 


As  v-»«,  the  interaction  coefficient 
C •  (x) -*0  and  expression  (5.8)  reduces 
to  the  exciting  force  obtained  in  the 
strip  theory  of  Salvesen,  Tuck  and 
Faltinsen  (1970).  The  short-wavelength 
limit  of  (5.7),  however,  does  not  lead 
to  the  same  result  since  the  slender- 
body  approximations  involved  in  (5.7) 
and  (5.8)  are  different. 

The  relative  error  involved  in  the 
inner  diffraction  problem  :s  a  factor 
1+0 (e2vcos6) .  The  corresponding  error 
factor  for  the  radiation  problem  is 
1+0 (c 7  )  for  all  wavelengths.  It  is  thus 
expected  that,  unlike  the  radiation 
velocity  potential  and  the  exciting 
force  (5,8),  the  accuracy  of  the  diffrac¬ 
tion  potential  and  the  corresponding 
exciting  force  (5.7)  decreases  with 
decreasing  wavelength.  This  is  support¬ 
ed  by  the  comparison  of  the  exciting 
force  and  moment  predicted  by  (5.7)  and 
(5.8)  with  an  exact  three-dimensional 
theory. 

The  equations  of  motion  follow  by 
equating  the  inertia  forces  to  the  sum 
of  the  pressure  forces  and  the  ship 
weight.  Combining  (5.1)  and  (5.4)  and 
taking  the  origin  of  the  coordinates 
on  the  ship  centerline  and  above  the 
center  of  gravity,  the  equations  of  the 
heave  and  pitch  motions  follow  in  the 
form 


j  —  3 , 5 


n -t. .  +  a  x 
3  13  i 


l5-9) 


where  M^pV  is  the  ship  mass,  M55  is 

the  ship  moment  of  inertia  with  respect 
to  the  y-axis  and  M35=M53=0  ^ue  to  the 

special  choice  of  the  coordinate  system. 
Substituting  (5.2)  in  (5.9)  and  rear¬ 
ranging  terms  we  get 

l  nj{-wJ(Mij+aij)+iwbij+C.j) 

3“  J  ,  3 


=  AXi, 


i=3, 5 . 


(5.10) 


This  is  a  system  of  linear  equations 
that  can  be  easily  solved  for  the 
complex  amplitudes  n..,  j=3,5. 

Numerical  computations  of  the  hydro- 
dynamic  forces  and  motions  were  per¬ 
formed  for  a  Series  60  hull  (C„=0.7, 
parent  form).  The  Salvesen,  et.  al. 

(1970)  strip-theory  results  are  also 
shown  together  with  an  exact  three- 
dimensional  numerical  solution  by  Inglis 
(1980)  and  experimental  data,  where 
available. 

The  results  for  the  radiation  prob¬ 
lem  are  shown  in  Figures  2  and  3,  and 
are  compared  to  experimental  data  of 
Gerritsma  (1966) .  The  agreement  be¬ 
tween  the  unified  theory  and  experi¬ 
ments  is  very  good  both  for  the  hydro- 
dynamic  force  distribution  (Figure  2) 
and  the  added-mass  and  damping 
coefficients  (Figure  3  adapted  from 
Newman  and  Sclavounos  (1980)),  and 
indicates  a  notable  improvement  over 
st'-ip  theory.  The  deviation  of  a5c 
from  the  experiments  that  occurs  at  low 
frequencies  is  supported  by  the  excel¬ 
lent  agreement  between  the  unified 
theory  and  the  exact  three-dimensional 
solution  of  Inglis  (1980)  for  all 
frequencies. 

The  slowly  varying  sectional  hydro- 
dynamic  force  distribution  for  the 
diffraction  problem  is  given  by 

F^ (x)=  -  iwpAC?  (x)  |  <)>7Hn3dl.  (5.11) 


The  corresponding  expression  for  the 
Salvesen  et  al.  (1970)  strip  theory  is 
obtained  from  (5.8)  with  C^(x)=0,  in 
the  form 


F3(x)  =  -^£3*  (iuH3-*3flj). 


evzcos (vysinB) dl ,  (5.12) 

where  all  the  quantities  in  (5.11)  and 
(5.12)  are  defined  in  Section  2B.  The 
results  for  Fi  and  the  amplitudes  of 
the  exciting  force  and  moment  shown  in 
Figures  4  and  5  respectively,  indicate 
that  the  strip-theory  predictions  are 
in  general  higher  than  those  of  the 


unified  theory. 

The  two  expressions  of  the  exciting 
force  and  moment  (5.7)  and  (5.8),  de¬ 
rived  respectively  by  pressure  integra¬ 
tion  and  using  the  Haskind  relations, 
compare  well  with  each  other,  being  in 
closer  agreement  for  smaller  values  of 
v|cosb|.  For  head  waves,  both  agree 
well  with  the  experimental  data  of  Vugts 
(1971)  and  the  three-dimensional  solu¬ 
tion,  with  a  relatively  more  favorable 
agreement  of  (5.8)  with  the  theory  of 
Xnglis  (1980) . 

Finally,  the  motions  of  the  ship  free 
to  heave  and  pitch  are  presented  in 
Figure  6,  where  comparisons  are  made 
with  three  different  sets  of  experi¬ 
mental  data  due  to  Nakamura  (19*6)  , 
Shintani  (1966)  and  Yamanouchi  and  Ando 
(1966).  For  the  model  tested,  the 
center  of  gravity  is  taken  at  the  origin 
of  the  coordinates  and  the  longitudinal 
radius  of  gyration  is  taken  equal  to 
1/4  of  the  ship's  length  L. 

All  theories  agree  well  in  general 
with  the  experiments,  except  for  short 
waves  at  6=180°.  The  conclusions  drawn 
for  the  exciting  forces  and  moments 
predicted  by  (5.7)  and  (5.8)  remain 
unchanged  for  the  respective  ship  motion 
amplitude  predictions,  since  the  added- 
mass  and  damping  coefficients  are  the 
same  in  the  two  versions. 


data.  The  strip  theory  predicts  well 
the  motion  amplitudes  and  phases  but  is 
in  less  satisfactory  agreement  with  the 
other  theories  and  experiments  for  the 
hydrodynamic  coefficients  and  the 
exciting  forces. 

The  integral  equations  (2.27)  and 
(2.34)  are  solved  by  iteration  and  the 
solution  obtained  in  this  manner  has 
been  checked  against  an  independent 
matrix  inversion  solution.  The  computa¬ 
tion  time  required  for  this  task  is 
minimal,  leaving  the  solution  of  the 
two-dimensional  problems  as  the  main 
computational  effort  involved  in  the 
unified  theory. 
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Figure  2  -  Longitudinal  distribution  of  the  added-mass  and  damping  sectional 
hydrodynamic  force  for  a  Series  60  hull  (Cg=0.7)  at  X/L=  1.706,  0.758  and  U=0 
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Figure  3  -  Heave  and  pitch  added-mass  and  dampinq  coefficients  for  a  Series  60 
hull  (Cjj=0 . 7 )  at  U=0  (  Adapted  from  Newman  and  Sclavounos  (1980) ). 
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Figure  4  -  Longitudinal  distribution  of  the  real  and  imaoinarv  part  of  the  sectional 
exciting  force  (with  the  exp(ivxcosB)  component  factored  out)  for  a  Series  60  hull 

(C  =0.7)  fixed  in  head  (6=180°)  and  bow  (6=135°)  waves  at  \/h=  1.706,  0.758  and  U=0. 
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